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Abstract 

The MC toy model for Bose-Einstein correlations (BEC) is considered to 
make a best choice from different reference distributions. It occurs that the 
minimal bias in the BEC parameters estimation is provided by the reference 
sample which is being simulated from "observed" sample by the turn of all 
Pt vectors in an event by a random 5^. 

It proposed to use similar approach in the analysis C2(At],A0) correla- 
tion. To make the analysis less dependent on reference sample, new defi- 
nition of A(/) is implemented as the angle between vectors p\ and p2 in the 
transverse plane. Then C\ (A0 ) correlation can be studied without a refer- 
ence sample subtraction. 

The ridge in C2(A77,A(/>) correlation appears if an appropriate cut in Q- 
value (2 = 7^2 — Pi) is apphed. 



1 Introduction 

One of the way to study the space-time picture of hadronic interaction is to mea- 
sure the correlations between two identical particles which momenta are close to 
each other. This correlations are used to estimate the size of the domain, from 
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which particles produced are emitted. 



Consider the amplitude M{ri,r2) in coordinate representation in which two 
identical pions, TTi and 712, are emitted at the points ri and r2 correspondingly. 
The coordinates ri, r2 can not be measured, and instead one does estimate pions 
momenta quantities pi and p2 ■ So we have to take the Fourier transform 

M„(p,.p.) = /(^.-'.^.-W(n..,). (1) 

Besides this we have to consider the permutation of two identical pions. That is 
we have to add to Ma the amplitude 

M»(p„/..) = /gl.-.g^.--M(.„..). (2) 

where the pion with momentum p2 was emitted from the point ri and wise versa. 
This can be written as 

M{pi,p2)=Ma+Mb = Ma-{\+e''Q), (3) 

where the 4- vectors r = ri — r2 and Q = P2 — Pi^- 



Finally the cross section takes the form 

d^pid^P2 2V "^^^ ^ \Ma\[i+<e >) 



(4) 



Here the factor 1/2! reflects the identity of two pions and the angular brackets 
indicates the averaging over the (ri, r2) space distribution. Assuming, for simplic- 
ity, the Gaussian form we get < e"'^ >= e^^'' . 

Thus the width of a peak at small momenta difference characterizes the 

size < r > of the domain from which the pions were emitted. 

The general definition of a correlation function of 2 variables - pi and p2 can 
be written as 

C2[Pl,P2) = T- (5) 

Pref{p\,Pl) 



Everywhere we suppose that particles ai^e pions 
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The value PrefiPii Pi) is called a reference function, artificial function with switched 
off the correlation under search . Usual task is to extract correlation C2 from the 
"measured" value P{pi,p2). In the following we will consider Bose-Einstein cor- 
relations (BEC) in the multi-particle final state (MPS) in pp-coUisions. Suppose 
that BEC correlation is determined by P{p\,p2) with the reference Prefipi^Pi)- 
Usually MC generators simulate MPS without BEC. In the following the ATLAS 
MC [T] of the pp-interaction with 7 TeV center of mass energy is used. In the MC 
event sample there is no BE correlations , but known hadronic resonance states 
are "produced". 

The distributions of identical bosons should be influenced by BEC. The prox- 
imity in phase space between final state particles with 4-momenta pi and p2 can 
be quantified by 
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Q = \l-{pi-P2V- (6) 

The BEC effect is observed as an enhancement at low Q ~ 300MeV . To extract 
the effect one can compare measured Q-spectra with similar one but without BEC, 
with reference spectra . Then the ratio 

clN _ dNi-ef 
dQ 

can be fitted with an appropriate formulae 

R(Q)=F(rQ)+a + bQ (8) 

F(rQ) is the Fourier transform of the spatial distribution of the emission region 
with an effective size (r). The function can be parameterize as a linear exponent 
F = Ce ''^^ or a Gaussian F = Ce^'s^ . The parameter b may be accounted for 
a long range correlation.For the same (r) there is quite simple relation between 
these two cases : 

re = y/nrg. (9) 

These two choices are connected to different radiation sources distribution inside 
the emission region. The exponential form approximates the case when radiation 
is originated from the surface of the sphere with radius r^. The second form 
describes the radiation in space ~ e^'^s^ . The fit of the same distribution 



3 



with these two cases gives r-parameter ratio as in ( |9). The exponential fit is a bit 
better and will be used in the following. 



2 MC BEC model and different reference functions 

Two values of ^ are measured in the experiment: for the like sign particle 
pairs and for the unlike sign case. To simulate BEC into MC distribution , 
each entry into can be weighted by the function 

W{Q) = \. + e-''Q (10) 

to produce ^^"^^ _ jj^g value r is fixed equal to 2fm. Of course, such model is 
the oversimplified toy but good enough to illustrate difference between reference 
functions. 

First consider the ideal case when MC will be used as the reference func- 
tion Fig [U shows the difference - and the result of the fit of the 
ratio dV]) with the formulae (|8]) [§. The result is somewhat trivial and can be consid- 
ered as a "goodness" of the MC generator. The procedure used might be looked 
as switching ON the BEC correlation, to add such correlation to MC procedure. 

In the BEC analyses of experimental data different approaches have been ex- 
plored and it is worth to compare known reference functions with the use of the 
toy MC model. 

The unlike Q-distribution looks as a natural choice because there should be no 
BEC for the unlike sign pairs. The source radius from the fit is close to the model 
value ( Fig [2l), however the resonances contribution makes fit quality bad. The 
contributions of the and the remnants of rj mesons are clearly seen. 

One has to find a way to switch OFF BEC from the "observed" ^^^^ distri- 

j\T ref 

butions. It looks as an easiest way to construct the Jq from two independently 
detected events: one momentum vector is taken from event under analysis, and 
the second - from a preceding events. The result is strongly dependent on require- 
ments of such event selection. Moreover, such requirements formulation is rather 
arbitrary. 



^ If the parameter a or is equal zero within two statistical errors it is fixed to be zero. The 
estimations of parameters error are "normalized" to make x^/Ndof = 1 
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Another approach is to emulate independent event sample with the use of the 
events set under consideration. The "referenced" event will have at least the same 
particle multiplicity,the same numbers of positive and negative charged particles. 
Evidently, such procedure will work as a blind correlation terminator because 
there is no any BEC "marker" on a particular entry. The hope can be that pairs 
with small Q-values will be replaced by pairs which had higher Q-value before 
the transformation. The reference distribution '^'^'i ^ will contain one vector from 
the real event and another from the reference one. 

The standard procedure to prepare mirror sample from the original one is to 
change in the event all momenta 3-vectors as p — )■ —p (PMIR case). It would be 
good if events (and the detector) have such symmetry. The reference sample is the 
result of measurements with the detector identical to the real one but in the mirror 
space. Fig [3] illustrate results of such procedure. We do see that the distribution 
'^^^Q ^ is different from the model too much. In particular there is no entries in the 
reference distribution if Q-value in the real events is smaller than ~ 300MeV : the 
procedure removes all small Q entries from the reference Q-distribution, not only 
BEC entries. In other words, the BEC contribution will be overestimated and the 
ratio dS]) Q-distribution will be wider ( compared with Fig [1} and the r-parameter 
of ([8]) will be smaller than in the model. 

The better way might be if the reflection is made only in the transverse plane 
Pt — )■ —pt (PTMIR case). The results are a bit closer to the model values but still 
there is no small Q events in the reference distribution if the model Q-values is 
smaller than ~ BOOMeV ( FigH). 

One gets a bit better reference distribution if vectors pt in an observed event 
will be tumed by a random value of : the same tum for all tracks - PHIRAND 
case , or each track is tumed by a random angle - PHIWIDE case ( Fig (6]). Fit 
parameters for the PHIWIDE case is found to be closer to the model than in all 
previous cases. 

The model we consider, the radiation zone radius has a constant value. How- 
ever one might suspect that "imperfection" of reference samples is larger at small 
particles multiplicity. This is the case indeed ( Fig |7]). The samples with the reflec- 
tion algorithms produce quite detectable multiplicity dependence. The approach 



'^The "measured" <p distribution is slightly non uniform because of real detector acceptance and 
detection efficiencies. The same is true for 7] distribution 
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with the random pt turn again looks better. 

Another important feature is the value of a systematic uncertainty. One can 
produce the reference sample of MC events without BEC and estimate a value 
of the source radius with different reference samples. It occurs that the value is 
quite small ( .15 fm, Fig ??) and might be consider as a contribution of non-BE 
correlation with small Q value or as a systematic uncertainty. In the following, we 
will see that non-BE correlation at small Q exists indeed. 

3 2- and 3-dimensional correlations 

Until now, we have used 1 -dimensional distribution - number of events as a func- 
tion of Lorentz invariant value Q. This is a natural approach for a simple fit to 
measure phenomenon parameters . Another way have to be considered in the case 
of a search of deviations from a phase space prediction with a smooth reaction 
amplitude. Usually 2-dimensional correlations is being used. Let us consider the 
traditional correlation plot for like sign pairs where differences 



A02;r = <|>2-01, ^^^^ 

^r\ = 7]2- r\i 

will be taken as arguments - C2{^(^2m^^)- Q- value and Nch ( number of 
charged particles in the event) can be used as additional arguments to form 3- 
dimensional correlations C3 (A02;r7 Arj , Q) and C3 (A^2;r7 At] , Nch) ■ 

Fig [8] shows several correlation plots. If particles would have the uniform 
distribution for 77(— 2.5,2.5) and for ^{—K,n) , then the values Ci(A7]) and 
Ci(A^2;r) will have triangle distribution in the range T](— 5.,5.) and ^{—In.ln). 
For these reasons C2(A02;i:, At]) looks as a pyramid. Except rather strong corre- 
lation at A^2n ~ Ar] ^ nothing significant can be seen . One have to subtract 
a reference distribution (one of the pt vector of the pair is turned by a random 
d(j) ) to see some additional structure. Unfortunately small nonuniformity might 
appear artificially because of imperfections of the reference sample. Some fea- 
tures are connected with definition of A^2;r: evidently,A^2;r = ±2;r is the same as 
A02;r = 0. 
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It is better to change the definition of 



A0;r = <^)2-<^>1,|<^>2-<^>1| < ^ ^^^^ 

A^jj: will have a uniform distribution in the range (—n, n) in the absence of cor- 
relations and to study correlations one has no need for any reference to subtract. 
Fig [9] shows C2(A0,Ar]) before and after reference subtraction. Similar plot for 
unlike sign pairs is shown also . Strong correlation is independent on charge of 
particles. 

Let us consider other correlations C2(A^,2) and €2(^7], Q) (Fig [TOl). An 
enhancement at small Q value can be expected. Q is small if A0 ~ and At] ~ ( 
momentum value of particles in pairs are larger than masses) ^independent on the 
momentum values and the charges of particles. At the same time at large Q value 
there is a dip at A0 ~ and ( Ar] ~ ). Then after integration in Q , there might be 
full compensation of the enhancement at small Q , if there is no dynamical reason 
to make an enforcement of this enhancement. By other words, with an appropriate 
cut in Q value, small A0 and At] correlation can be suppressed (Fig [TT}. An effect 
of dynamical correlations exists only when Q is less than ~ lOOMeV . After such 
cut nonuniformity in A^ decreases from ~ 4% to ~ 1%. Very important that such 
behavior is the same for like sign and unlike sign pairs ( except natural traces of 
resonances at 2 > 400 MeV) . 

The next step is to have a look an influence of such Q cuts on 2-dimensional 
plots (Fig [T2I) . The "Long-Range, Near-Side Correlations" (ridge) starts to be ap- 
pear. Such effect has not been seen in MC before ( see [i3 |)not only because Q-cut 
never was used before, but also because of different way to construct the reference 
sample ( see in the follows). The ridge was observed [3] in the data analysis of pp- 
coUisions at the LHC at 7TeV in high multiplicity events with selection of pairs 
of charged particles with pt in the range (1-3) GeV. In MC distribution there is no 
trace of a ridge, if only small pt cut has been done. Fig [13] shows the correlations 
for particles with pt > IGeV . The multiplicity cut only decrease statistics but does 
not change the conclusion. By other words, the selection of pairs of particles with 
Q > 200 MeV can be good tool for search of a ridge effect. 
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4 Conclusions 



Rather simple MC model for Bose-Einstein Correlations is used to judge the qual- 
ity of different ways to extract the BEC from observable distribution. Several emu- 
lated reference samples have been considered. The reference sample with random 
turn of observed transverse momentum vector pt randomization ) provides the 
fit results closest to the model values. 

Such approach is implemented for the analysis of C2(A0, At]) correlation. In 
the traditional definition of A02;i:, to extract correlations the subtraction of a ref- 
erence distribution is the must. New definition A^^r as the angle between pti and 
Pt2 makes A^;j distribution uniform with small contribution from an intrinsic cor- 
relations. For this reason the Ci(A^;i:) can be studied without any reference. It 
is shown that at small Q- value there might be two kind of Ci (A^j^) correlation - 
kinematical and dynamical. At large Q there is Ci (A^^i) anticorrelation. Then so 
called short-range correlation ( A^^r ~ ) can be suppressed by an appropriate cut 
in the Q- value. 
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Figure 1 : MC Model fit with the "ideal" reference 
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Figure 2: MC Model fit with the 4^ for unlike sign pairs as the reference 
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Figure 3: MC Model fit of the ratio R{Q) where the reference sample consists of 
real events with all 3-vector momentum reflected in space 
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Figure 4: MC Model fit of the ratio R{Q) where the reference sample consists of 
real events with all pt momentum reflected in the transverse plane 
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Figure 5: MC Model fit of the ratio R(Q) where the reference sample consists of 
real events with all pt momentum turned by the same random 50 
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(b) Fit of the ratio = 



Figure 6: MC Model fit of the ratio R(Q) where the reference sample consists of 
real events with all pt momentum turned by the same random 50 
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Figure 7: Fit results for different reference samples 
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Figure 8: MC C2(A(^, At]) for like sign pairs correlation, = A^27i 
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(b) MC 7?(A0,A77) = c^(A^.A,)-cf (AfA,) 
like sign pairs 



Figure 9: MC C2(A(^, A?]) ,A(^ = A^jt 
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(a) MCC2(g, A;]) correlation 
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(b) MC C2(2, A(^) unlike sign correlation 



Figure 10: MC C2(2, A^) and C2(2, At]) correlations 
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Figure 1 1 : MC Ci (A0) distribution with different cuts in Q value 
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(a) MC R{A(j),Ari) correlation without Q cut 




(b) MC R{A(j),A'n) correlation with Q > 120 MeV 
cut 




(c) MC R{A(l),Ar]) correlation with Q > 200 MeV 
cut 



Figure 12: MC i?(A^, At]) correlation 
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(a) MC R{A(j),Ari) correlation without Q cut 
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(b) MC 7?(A0,At7) correlation with Q > 120 MeV 
cut 
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(c) MC R{A(l),Ari) correlation with Q > 200 MeV 
cut 



Figure 13: MC i?(A(^, At/) correlation for pairs with p'^'" = IGeV 



21 



